Testing the topological nature of the fractional quantum Hall edge 
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We carry out numerical diagonalization for much larger systems than before by restricting the 
fractional quantum Hall (FQH) edge excitations to a basis that is exact for a short-range interaction 
and very accurate for the Coulomb interaction. This enables us to perform substantial tests of 
the predicted universality of the edge physics. Our results provide compelling evidence that the 
behavior of the FQH edge is intrinsically nonuniversal, even in the absence of edge reconstruction, 
and therefore does not bear a sharp and unique relation to the bulk FQH state. 
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PACS numbers: 

The interior of a fractional quantum Hall (FQH) sys- 
tem [1 is gapped, but massless excitations exist at its 
edge, which constitutes a reahzation of a one-dimensional 
electron liquid described generically by the Tomonaga- 
Luttinger theory [?. Much attention has been focused 
on the edge physics since the work of Wen [3^, where 
it was conjectured that the exponent describing the 
long distance, low energy physics of the chiral (unidi- 
rectional) FQH edge is a unique "topological" quan- 
tum number for any given FQH state, independent of 
details, just as the quantum Hall resistance. Our un- 
derstanding of the ordinary one-dimensional liquids is 
largely based on the method of bosonization, which ex- 
ploits a one-to-one correspondence between the fermionic 
and bosonic Fock spaces in one dimension, and identi- 
fies a relationship between the operators on these spaces; 
specifically, the fermionic field operator 'ip{x) is related 
to the bosonic field operator 4>{^) through the expres- 
sion '4^{x) ~ exp[— i(/)(a:)], which can be established rig- 
orously at the operator level [5]. In the absence of a 
similar rigorous derivation for the electron field operator 
at the edge of a FQH system, Wen formulated an effec- 
tive field theory approach (EFTA) fT wherein he postu- 
lated that the electron operator at the edge of the 1/m 
FQH state, defined by Hall resistance quantization at 
Rh = h/{l/m)e^, is given by 
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Antisymmetry under exchange quantizes m to an odd 
integer value, independent of parameters other than the 
quantized Hall resistance, which leads to universal prop- 
erties for the edge physics. A direct test of this asser- 
tion is through tunneling of an external electron later- 
ally into the edge of FQH system. For the fractions 
V = n/{2np -I- 1), a generalization of Eq. (IT]) predicts 
the I-V characteristic for electron tunneling into a Fermi 
liquid to be / '^ V^ , where the tunneling exponent has a 
universal value of a = 2p + I. Ingenious experiments 
0HZ] have measured the edge exponent by determin- 
ing the I-V characteristics for tunneling from a three- 
dimensional Fermi liquid into the FQH edge. While they 
establish the existence of non-Fermi liquid (Tomonaga- 



Luttinger) behavior, with an exponent different from a 
one-dimensional Fermi liquid (a = 1), they also show 
discrepancy from the EFTA prediction of Wen. Specifi- 
cally, experiments find an exponent that is continuously 
varying with the filling factor, and thus is not determined 
solely by the quantized Hall conductance. Furthermore, 
the measured exponents are significantly different from 
the EFTA predictions: at filling factors ly = 1/3, 2/5, 
and 3/7, the exponents are ^ 2.7, 2.3, and 2.1, respec- 
tively [4~7i], to be compared to the EFTA prediction of 
3.0. 

A number of theoretical papers have addressed this 
inconsistency [7H20] . Some of these suggest that the dis- 
agreement is due to edge reconstruction, which produces 
several counter-propagating edge modes (for which the 
exponent is not universal) P"51[T51 ^T\ , while some pro- 
pose that the inconsistency persists even in the absence 
of edge reconstruction, thus pointing to a more funda- 
mental deficiency of the EFTA |13l [19] . A resolution of 
this issue is important not only in its own right, but also 
in view of the potentially useful notion that the character 
of a bulk FQH state can be ascertained from the behav- 
ior of its edge physics 22 , which logically rests on the 
existence of a unique relationship between the two and 
hence the universality of the latter. 

Exact diagonalization studies often provide an unprej- 
udiced, reliable, and decisive tool for testing ideas in the 
field of the FQH effect. For the edge physics, however, 
it has not been clear if the discrepancy between the fi- 
nite system results and the EFTA is intrinsic or a finite 
size artifact; finite size corrections are more severe for 
the edge physics [20] because of power law decay of cor- 
relations, in contrast to the Gaussian decay in the bulk. 
Unfortunately, the dimension of the Hilbert space grows 
exponentially with the number of electrons, making it 
impossible to increase the system sizes significantly in 
exact diagonalization studies. 

In this Letter we report on substantial microscopic 
tests of the EFTA by diagonalizing the Coulomb Hamil- 
tonian in a truncated space of edge excitations. Specif- 
ically, we consider the edge excitations of the 1/3 FQH 
state in the disk geometry, and the truncated space con- 
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FIG. 1: [Color online] Energy spectrum for the edge exci- 
tations of u = 1/3 for A'' = 9, 27, 45 particles at electron- 
background separations in the range d = 0.0 and 2.5. Blue 
dots indicate the energies obtained by CF diagonalization, 
whereas the adjacent red triangles (shifted along the x axis 
for clarity) show the bosonic spectra (see text for explana- 
tion). All energies are quoted in units of e'^/eZ, and measured 
relative to the energy of the ground state at AM — 0. AM 
is the angular momentum of the excited state. 



tains all states of the form 

'^a=l[('^i-^^kf^f, Ad = M*+NiN-l), (2) 

j<k 

where Zj = Xj — iyj denotes the electron coordinates as 
a complex number, M is the total angular momentum 
of the sate, and <I>^ are all lowest Landau level states 
(labeled by a — 1, • • • , D*) at total angular momentum 
M* . The dimension of this basis space is much smaller 
than the full dimension of the lowest Landau level states 
at M, which makes it possible to investigate much larger 
systems; we have studied as many as 45 particles. The 
restriction to this basis is equivalent to restricting com- 
posite fermions [23] to their lowest A level (also known as 
composite fermion Landau level) ; the space of states can 
be enlarged in the standard manner [21] by also includ- 
ing at M* states occupying successively higher A levels 
(and projecting the total wave function onto the low- 
est Landau level), but that will not be necessary for our 
present purposes. The lowest- A-level approximation for 
composite fermions is known to be excellent, and we have 
also confirmed its accuracy explicitly for edge excitations 
for systems with six and seven particles, for which exact 
results are available, both with and without the confine- 



ment potential. Also, Vf^ are the only states that sur- 
vive if we add to the Coulomb interaction an appropriate 
infinitely strong short range interaction that annihilates 
states containing electronic pairs with angular momenta 
equal to unity; our results below are exact for this model. 
Therefore, we believe that our truncated Hilbert space 
ought to capture the topological nature, if it exists, of 
the edge physics. We believe that this model actually 
gives the best chance for universal behavior; mixing with 
higher A levels can only spoil it [T5] . 

We consider a system of two dimensional electron gas 
in disk geometry. The neutralizing background has uni- 
formly distributed positive charge contained in a disk Q^ 
of radius Rn = \/2N/v for a system of N particles at 
filling factor v; the positively charged disk is separated 
by a distance d from the electron disk (quoted in units of 
the magnetic length below). The electrons are approx- 
imately confined to the same radius because of charge 
neutrality in the interior. This system is modeled by the 
following realistic Hamiltonian: 



H = Ek + V,, + Kb + H 



bb 



(3) 



where the terms on the right hand side represent 
the kinetic, electron-electron, electron-background, and 
background-background Coulomb interaction energies, 
respectively. At large magnetic fields only the lowest 
Landau level states are occupied, hence the kinetic energy 
Tiujc/2 (where Wc = eB/nibC is the cylcotron frequency) 
is a constant and will not be considered explicitly. 

The wave functions 4"^ are in general not orthogo- 
nal, and we use the method of composite-fermion di- 
agonalization (CFD) [25 to orthogonalize them by the 
Gram-Schmidt procedure, evaluate the Hamiltonian ma- 
trix elements, and diagonalize it to obtain the eigenval- 
ues and eigenvectors. All matrix elements and scaler 
products needed for this purpose are evaluated by the 
Monte Carlo method, as explained elsewhere in the lit- 
erature [24j. While sufficiently accurate energy spec- 
trum requires ~10-20 million Monte Carlo iterations, 
the spectral weights require ^ 200 million iterations for 
each eigenstate. These numbers do not vary significantly 
with N, but the computation time increases exponen- 
tially with N and AM, limiting our study to systems 
with A^ = 45 for energy, and to A^ = 27 for spectral 
weights. The energies were calculated for AM = 1 — 8 
and the spectral weights for AAf = 1 — 4. 

Using the CF diagonalization procedure, we compute 
the spectra of edge excitations of the 1/3 state, shown 
in Fig. (II]), for several parameters in the range A^ = 6- 
45, d = 0-2.5, and AM = 0-8. Our large system cal- 
culations confirm an earlier study [TS] that edge recon- 
struction occurs for d larger than a critical separation, 
which is approximately 1.5-2.0 magnetic lengths. For d 
greater than the critical separation, a simple explana- 
tion for the observed nonuniversality follows in terms of 
edge reconstruction; a model that assumes a single chiral 



mode is inadequate to describe experiments, forcing one 
to consider multiple edge modes, which produces nonuni- 
versal results. However, the important question remains 
whether universality occurs in the absence edge recon- 
struction. 

We address this issue by following the pioneering work 
of Palacios and MacDonald,^ to test the validity of Eq. 
M, upon which the notion of universality rests. Specifi- 
cally, we compare certain matrix elements of the electron 
field operator, computed from our CF diagonalization re- 
sults, with the predictions of the bosonized form in Eq. 
(1). We also consider d greater than the critical separa- 
tion for completeness; here, we assume that the ground 
state remains at AM = 0, which can be arranged by 
adding an ad hoc angular momentum dependent single 
particle energy term that strongly penalizes the edge ex- 
citations responsible for edge reconstruction, but does 
not change either the eigenfunctions or the energy or- 
dering of states at a given AM. (This can be accom- 
plished by adding an appropriate parabolic confinement 
term which adds to the total energy a term proportional 
to the total angular momentum.) 

The spectral weights are defined by 
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(4) 



where |{?T.i}) represents the bosonic state with occupa- 
tion {ni}, |0) is the vacuum state with zero bosons, 
and V'^(^) is the electron creation operator at posi- 
tion 9 along the edge circle. Here I denotes single 
boson angular momentum; the total angular momen- 
tum is denoted by AM — J^i ' "■( ^-^d the total en- 
ergy by AE = X); "-i^ii with e/ being the energy of 
a single boson at angular momentum I. With the 
help of ip^O) ex e' 
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4>-{0)^ it is straightfor- 



ward to obtain the predictions for the spectral weights: 
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We note that the denominator in Eq. Q eliminates the 
unknown normalization constant ^/z. 

In order to obtain the spectral weights from our elec- 
tronic spectra, it is natural to identify the vacuum state 
|0) with the ground state of interacting electrons at 
1/ = 1/m, denoted by |^^), and the field operator has 
the standard meaning oi'l|'^9) = J2i ''li{(^)'^i — J2i '/'i (^)i 
where c] and ci are creation and annihilation operators 
for an electron in the angular momentum I state. The 
denominator of Eq. Q corresponds to 
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FIG. 2: [Color online] A'^ dependence of spectral weights 
|C{„j}p for several states {n;}, indicated on each panel, and 
various separations d (quoted in units of the magnetic length 
in the {0100} panel). The EFTA prediction from Eq. (1) is 
indicated by a star on the y-axis, with the value also given 
on each panel. The points on the y-axis are determined by a 
quadratic fit to the finite A'' results. 



where |^^) is the ground state of N interacting electrons 
at V = 1/m, and Lq = mN . The numerator is similarly 
defined as as 



{{ni}\^He)\Q) = 



Here, we have 
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(6) where A is the antisymmetrization operator, A/l is the 



normalization constant, and L = Lq + AM is the angular 
momentum of added electron. 

The wave function ^r„^i^, the electronic counterpart 
of the bosonic state Kn;}), clearly represents an excited 
state at total angular momentum M — AM + mN{N + 
l)/2, which should also be related to the total angu- 
lar momentum of the N particle ground state through 
M = L+mN{N-l)/2. At each AM, there are in general 
many eigenstates. Following Ref. [TS] we identify ci, the 
energy of a single boson with angular momentum I, with 
the lowest energy at I = AM in the calculated spectrum. 
Using the equations J2i ^"-i = AAf and -E{n,} = J2i '^i^i^ 
the energies of the all bosonic states {ni} can now be 
obtained (see Fig. IT]), which can then be identified with 
the corresponding electronic states. We note that in the 
lowest A level subspace, the numbers of electronic and 
bosonic states are equal at each AM, so a one to one 
correspondence between the two sets of states can be es- 
tablished from their energy orderin. For small systems 
(for example, iV = 9 in Fig. [l]), the CFD spectra and the 
bosonic spectra are very close to each other, which ex- 
plicitly confirms the interpretation of the lowest branch 
as the single boson branch. The agreement between the 
electronic and bosonic spectra becomes less accurate with 
increasing N or AM, but still remains adequate for the 
low energy states, which will be our focus. (The higher 
energy states of the spectra shown in Fig. 1 mix with 
higher A level excitations of composite fermions, not con- 
sidered in our model.) 

Figure ^ shows the squared spectral weights for dif- 
ferent excited states as a function of N and d. A 
quadratic fit extrapolates the result to the thermody- 
namic limit 1/A^ — 0. The EFTA predictions from 
Eq. (1) are also shown in each panel. These plots demon- 
strate the central result of our work: the spectral weights 
are nonuniversal; they depend on d] and they do not ex- 
trapolate to the EFTA value. For the {1000} excitation 
the thermodynamic result agrees with predicted result of 
3.0 for all d; however, in this case our truncated space 
contains a single (center-of-mass) excitation, the wave 
function for which is independent of interactions (within 
our model), and therefore the agreement is not meaning- 
ful. For many cases, the deviation from the EFTA value 
is substantial; even the spectral weights of single boson 
states, such as {0100} and {0010} exhibit significant d 
dependence. 

Our study thus indicates that Eq. 1 is not valid for the 
1/m FQH edge for the Coulomb interaction, and there- 
fore there is no reason to expect the edge exponent to be 
a topological quantum number for the \/m FQH state; 
this conclusion very likely holds for other FQH states as 
well, given that their edges are believed to be more com- 
plex. The problem of how in reality the electron field 
is related to the bosonic field remains unresolved, how- 
ever. Following Ref. [12] one may abandon the anti- 
symmetry requirement and try an expression of the type 



"ipix) ~ g-^^/a(t>{x) -^j^jj arbitrary a; we have found that 
no single value of a gives a satisfactory description of 
all spectral weights that we obtain from numerical di- 
agonalization. It is also worthwhile here to mention the 
possibility of a nonlocal relation between the two, as sug- 
gested in Ref. [13] . 

Other models for edge confinement have been used. 
One such model [15] restricts the single particle angular 
momentum to a maximum value of Imax = 3(A^ —1)-\-Iq, 
which may be a reasonable approximation for cleaved 
edge overgrowth [7j. However, angular momentum con- 
servation shows (and our Monte Carlo calculations ex- 
plicitly confirm) that in this model the actual spectral 
functions identically vanish for AM > Iq, resulting in 
an even more substantial disagreement with the EFTA 
predictions. 

In summary, our study allows an estimate of certain 
thermodynamic properties of the FQH edge for a real- 
istic Coulomb model, and makes what we believe to be 
a compelling case that the behavior at the FQH edge is 
intrinsically nontopological, as also demonstrated by ex- 
periments. This result has obvious implications for the 
program of determining the nature of a bulk FQH state 
by probing its edge. 
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